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Wave  Propagation  in  Particulate  Media 


The  purpose  of  this  project  was  to  study  optical  pulse  propagation  in  a 
medium  consisting  of  random,  discrete  scattering  centers  with  a  view  to 
establishing  the  effects  of  clouds  and  fog  on  optical  communication  systems. 
The  delay  in  the  arrival  time,  duration,  and  amplitude  of  pulses  transmitted 
through  the  medium  and  the  spatial  spreading  transverse  to  the  propagation 
direction  were  to  be  measured  and  compared  to  the  predictions  of  theoretical 
siodels. 

An  appropriate  scattering  medium  for  experiments  of  this  type  was  found 
to  be  emulsions  of  two  dielectric  liquids  since  the  indices  of  refraction  of 
both  the  scatterers  and  the  bulk  medium  could  be  varied  by  choosing  different 
substances  and  the  size  of  the  droplets  of  the  dispersed  phase  could  be 
controlled  by  the  method  of  preparation.  In  addition,  the  emulsions  could  be 
stabilized  with  surfactants  and  have  much  higher  number  densities  than  can  be 
maintained  in  cloud  chambers.  This  allowed  simulation  of  the  effects  of 
kilometers  of  cloud  in  a  laboratory  scale  scattering  experiment. 

The  scattering  media  used  for  this  project  were  emulsions  of  either 
mineral  oil  or  diiodomethane  (DIM)  in  water  stabilized  with  lauryl  sulphate  as 
a  surfactant.  The  latter  emulsion  is  of  particular  interest  since  the 
relative  index  of  refraction,  1.74/1.33  *  1.30,  is  nearly  that  of  water  in  air 
and  thus  simulates  scattering  in  clouds  and  fog.  The  preparation  technique 
involved  forming  a  coarse  emulsion,  passing  it  through  nuclepore  filters  to 
refine  the  droplet  size  distribution  and  finally  exploiting  the  differential 
drift  velocity  under  the  influence  of  gravity  of  droplets  of  different  radii 
to  further  reduce  the  variance  of  the  size  distribution.  It  was  found 
practical  by  this  means  to  produce  tens  of  liters  of  the  scattering  medium 
with  well  controlled  parameters. 

The  parameters  determining  the  single  scattering  characteristics  of  these 
emulsions  were  determined  from  direct  measurement  of  the  index  of  refraction 
of  the  liquids  and  the  size  distribution  of  the  droplets.  The  diameter  of  the 
droplets  was  obtained  by  drawing  samples  from  the  scattering  media,  taking 
microphotographs  at  400X  and  measuring  the  size  of  the  images.  Since  the 
volume  of  the  sample  viewed  under  the  microscope  was  known,  the  number  density 
of  the  scatterers  and  the  complete  size  distribution  could  be  constructed. 
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The  mean  and  standard  deviation  of  the  droplet  diameters  of  the  systems 
studied  ranged  from  d  *  1.86  Pm,  ■  0.85  la  to  d  ■  12.9  Pm,  ■  3.8  Pm. 

The  number  density  of  the  scatterers  ranged  from  3.0  *  10 12  to  1.1  x  101**  m  3. 

The  mean  scattering  cross  section,  the  mean  of  the  cosine  of  the 
scattering  angle  and  the  average  single  scatter  albedo  were  calculated  using  a 
computer  code  for  Mie  scattering  for  each  of  the  prepared  emulsions.  The 
albedo  for  all  the  oil/water  systems  was  unity  to  six  significant  figures 
while  for  the  DIM/water  system,  the  lowest  calculated  albedo  was  0.99926.  The 
inverse  of  the  mean  scattering  cross  section  times  the  number  density  is  the 
scattering  length  or  optical  mean  free  path,  b  *  (Nog)  1 .  For  the  systems 
studied,  b  ranged  from  0.45  to  3.40  ran.  The  asymmetry  factor  g  *  1  -  <cos  0> 
ranged  from  0.050  to  0.171.  The  range  of  cross  sections  and  number  densities 
available  enabled  measurements  to  be  made  on  systems  trtiose  optical  thickness 
varied  from  0.1  to  170  all  in  a  scattering  cell  12  cm  long  and  40  cm  in 
diameter. 

The  multiple  scattering  experiments  performed  involved  the  measurement, 
with  a  streak  camera,  of  the  time  delay,  stretching  and  attenuation  of  532  nm 
laser  pulses  initially  of  25  ps  duration.  The  pulses  were  produced  by 
frequency  doubling  the  output  of  a  passively  modelocked  NdrYAG  laser.  Figure 
1  shows  a  typical  streak  camera  record  for  a  pulse  transmitted  through  4  cm 
(optical  thickness  t  *  63)  of  DIM/water  emulsion. 

Figure  2  shows  a  plot  of  time  delay  as  a  function  of  optical  thickness. 
The  experimental  data  represents  measurements  made  on  two  DIM/water  emulsions 
and  one  oil/water  system.  The  straight  line  is  the  behavior  predicted  by 
Monte  Carlo  simulation  calculations.^  The  DIM  system  data  agrees  much 
better  with  the  theory  than  does  that  of  the  oil  system  which  has  much 
stronger  forward  scattering,  g  ■  0.078  compared  to  g  *  0.142  and  0.171  for  the 
DIM  systems,  indicating  that  the  theory  and  simulation  results  do  not 
adequately  account  for  extreme  forward  scattering. 

Figure  3  displays  the  ratio  of  the  integrated  intensity  of  the  scattered 
pulses  to  that  of  the  incident  pulse  for  four  DIM  systems.  The  dashed  line  is 
the  Beer-Lambert  exp(-x)  and  the  solid  line  the  computer  simulation  result, ^ 
1.69  (gt  ♦  1.42)  multiplied  by  10  .  This  factor  of  10  discrepancy 

between  the  measured  and  calculated  values  is  unresolved  but  may  be  explained 
by  the  fact  that  the  detector  was  20  cm  from  the  surface  of  the  scattering 
cell. 


2 


0 


250  500  750  1000  1250  1500  1750 


Figure  1 


TIME  (ps) 


Typical  intensity  versus  time  streak  camera  record.  The  pulse 
on  the  left  is  a  reference  pulse  which  has  bypassed  the  scattering 
cell.  The  pulse  on  the  right  has  propagated  through  63  optical 
thicknesses  of  a  diiodomethane/water  emulsion. 


NORMALIZED  DELAY 


Figure  2.  Delay  in  the  mean  arrival  time  due  to  multiple  scattering. 

Measured  delay  times  have  been  normalized  by  multiplying  by 
c  b”1  g"°*94/0.62,  where  c  is  the  speed  of  light,  b  the 
scattering  length  and  g  the  asymmetry  factor.  °,A  -  DIM 
emulsions;  D  oil  emulsion.  The  line  is  t1*94. 


Relative  Integrated  intensity  versus  optical  thickness  times  the  asymmetry 
factor,  gT,  for  four  DIM  emulsions. 


The  general  conclusion  which  can  be  drawn  from  this  work  is  that  the 
computer  simulation  models^  of  multiple  scattering  of  optical  pulses  is 
adequate  for  describing  propagation  through  clouds  and  fog  and  other  systems 
with  moderate  forward  scattering  g  >  0.1.  Except  perhaps  that  the  actual 
energy  loss  may  be  somewhat  greater  than  predicted  by  the  model.  Systems 
exhibiting  extreme  forward  scattering  g  <  0.1  are  not  well  described  by  the 
theoretical  model.  A  more  sophisticated  model  which  takes  into  account  higher 
moments  of  the  scattering  function,1 2 3 4 5  e.g.,  <cos  0>  might  give  better 
agreement  with  the  experiments.  However  most  Mie  scattering  codes  currently 
in  use  including  the  one  available  for  this  program  do  not  evaluate  higher 
moments  making  comparisons  with  the  PSR  model!  difficult. 

A  complete  description  of  the  project  and  all  results  including  some  beam 
spread  and  field  of  view  effect  measurements  are  contained  in  a  technical 
report  to  be  published  in  Applied  Optics  under  the  title  "Multiple  Scattering 
of  Optical  Pulses  in  Scale  Model  Clouds."  Parts  of  this  work  were  also 
reported  at  the  Annual  Meetings  of  the  Optical  Society  of  America  at 
Kissimmee,  Florida  October  19,  1981-3  and  Tucson,  Arizona,  October  21, 

,  1982 at  the  Conference  on  Lasers  and  Electro-Optics,  Phoenix,  Arizona, 

April  16,  1982; ^  and  at  the  OSA  Topical  Meeting  on  Optical  Techniques  for 
Remote  Probing  of  the  Atmosphere,  Incline  Village,  Nevada,  January  11, 
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